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The flux tube between a quark and an antiquark in Coulomb gauge is imagined in the gluon- 
chain model as a sequence of constituent gluons bound together by Coulombic nearest-neighbor 
interactions. We diagonalize the transfer matrix in SU(2) lattice gauge theory in a finite basis of 
states containing a static quark-antiquark pair together with zero, one, and two gluons in Coulomb 
gauge. We show that while the string tension of the color-Coulomb potential (obtained from the 
zero-gluon to zero-gluon element of the transfer matrix) overshoots the true asymptotic string 
tension by a factor of about three, the inclusion of a few states with constituent gluons reduces 
the discrepancy considerably. The minimal energy eigenstate of the transfer matrix in the zero-, 
one-, and two-gluon basis exhibits a linearly rising potential with the string tension only about 
1.4 times larger than the asymptotic one. 
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1. Color- Coulomb potential and gluon-chain model 

The color-Coulomb potential, i.e. the ^-dependent part of the energy of a physical state of 
a heavy quark-antiquark pair at distance R in Coulomb gauge, was shown to represent an upper 
bound on the true static quark-antiquark potential [pj]. Numerical simulations demonstrated that it 
was asymptotically linear and its string tension, Gcoui> was measured to be 2-3 times larger than 
the standard asymptotic string tension, (7 asymp 0,0]. 

In this context, a series of questions naturally arises: How do flux tubes form in the Coulomb 
gauge? What mechanism is behind the collimation of color-electric fields into flux tubes? How 
does it reduce <7coui to a asymp ? 

In the gluon-chain model, proposed first by Tiktopoulos [Q] and later developed by Thorn and 
one of the present authors (J.G.) [a], the flux tube between a quark and an antiquark is visualized 
as a sequence of gluons. As a heavy quark and antiquark move away from each other, a chain of 
constituent gluons arises between them. The constituent gluons are bound together by Coulombic 
nearest-neighbor interactions. Schematically (see Fig. [I]): 

l^lf n > = 9(0) {ao + aiA + a 2 AA + a^AAA + . . . } q{R) |¥ ). (1.1) 



The goal of the present study was to test, in first-principle numerical lattice Monte Carlo 
simulations, the conjecture that constituent gluons do reduce the magnitude of the static quark- 
antiquark potential, and to measure the constituent gluon content of the QCD flux tube in Coulomb 
gauge. (For an early attempt in the same direction see Ref. ]^].) We will mainly report results of 
simulations of SU(2) lattice gauge theory at j8 = 2.4 on a 22 4 lattice. A complete set of our data is 
contained in Ref. fflb. 
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Figure 1: The gluon-chain model. 
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2. Transfer matrix and static potential 

The (Euclidean) time evolution in lattice gauge theory in Coulomb (or some other physical) 
gauge is governed by the transfer matrix 

= exp[-Jfa], (2.1) 
where is the Hamiltonian in that gauge, and a the lattice spacing. We used its rescaled form: 

T =(d^= expH ^- £ » R {12) 

where *¥q denotes the ground state and E$ is its energy. 

In an ideal case one would like to diagonalize T in the (infinite-dimensional) subspace of states 
containing a static quark-antiquark pair separated by a distance R. Then, the static quark-antiquark 
potential (in lattice units) could be computed from 

V(R) = -log(w), (2-3) 

where T max is the largest eigenvalue of the rescaled transfer matrix T in this subspace, correspond- 
ing to the minimal-energy eigenstate. In reality, this program cannot be realized; instead one has to 
reduce the subspace of states considered to a manageable size. Fortunately, if the quark and anti- 
quark are not too far apart, one can expect amplitudes of states with a large number of constituent 
gluons to be negligible, and can seek for the minimal energy eigenstate in a sector containing only 
the quark-antiquark pair plus a small number of constituent gluons. 

So we will diagonalize T in a finite M-dimensional subspace of trial "chain" states: 

\V k )=q a (0)Q?[U]q b (R)\ x ¥o) k=l,2,...,M, (2.4) 

where Qu are gluonic operators, functionals of the lattice gauge field, that depend on quark/anti- 
quark positions and some number of variational parameters (see the next section for a particular 
choice of the operator basis); a and b are color indices. All quantities of interest can be estimated 
from matrix elements 



f*, = <WF/> = (^Tr[fiJ(/ + 1)^(0,06/(0^0(^0 



On = OPtHF,) = (±Tr [gJ(0Gz(0] ) - ( 2 -5) 

■ , (2 " 6) 

computed by lattice Monte Carlo simulations. 1 Knowing these matrix elements, one can construct, 
via the usual Gram-Schmidt procedure, an orthonormal set of states = 1,2, . . . ,M}, then 

the matrix elements 

T kl = (<S> k \%\<S>i), (2.7) 

and finally determine the largest eigenvalue T max of the M x M matrix T. Such a calculation 
has to be repeated for various variational-parameter sets to determine the one which minimizes 
— log(T max ) at a given R. An estimate of the static potential in the sector of variational gluon-chain 
states, nicknamed below the "gluon-chain potential", will then be given by 

Vchain(tf) = -log(w)- (2-8) 



'The notation Q(t) indicates that the operator Q[U] is evaluated using links on a hypersurface of fixed time t. 
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3. Choice of the operator basis 



Our choice of the operator basis in Eq. ( [2.4] ) was dictated mainly by simplicity, and by some 
amount of trial and error. 

A gluon chain is assumed to consist of a certain number of constituent gluons between the 
heavy quark and antiquark, and the gluon ordering in color indices is correlated with their spatial 
positions between the heavy color sources. As usual, in a variational approach the optimal energy 
states represent a compromise between kinetic energy and interaction energy. While the kinetic- 
energy contribution prefers spatial derealization of gluons, the Coulombic interaction energy fa- 
vors as small as possible transverse displacement from the line connecting the quark and antiquark 
positions. To satisfy both requirements, the derealization in the spatial direction was achieved by 
a superposition of gluon operators in along the line joining the sources, while derealization 
in transverse directions could be realized by using "transverse-smoothed" gauge-field operators on 
the lattice, in which high-frequency components of the field are (e.g.) Gaussian-suppressed in the 
directions transverse to the qq-lme. The transverse smoothing introduced a single parameter p , the 
only variational parameter used in our operator Ansatz (see below). To further simplify this pilot 
study, we restricted the number of constituent gluons to at most two. 

Our procedure was thus the following (for further details see Ref. 

• The lattice configurations were fixed to the Coulomb gauge by standard methods. From 
lattice link matrices U we constructed 

Ai(x,t) = ±[Ui(x,t)-U?(x,t)\, B i (x,t) = l-±Tr[U i (x,t)}. (3.1) 

• These quantities were then Fourier-transformed, and we suppressed high-momentum com- 
ponents in directions transverse to the line joining the qq pair (e.g. the j-th direction): 

A f (k,t) ^exp [-p(k 2 -k 2 j)] A,(k,f), ->exp [-p(k 2 -^ 2 )] Bi (k,t), (3.2) 

then transformed them back to coordinate space, to get A- (x,f) and B/ (x,f), the A- and B- 
fields smeared in directions transverse to tj. p is a variational parameter, used to maximize 
the largest eigenvalue of the transfer matrix in a chosen basis of states. 



It was also useful to define, for i ^ j, the averages: 



AP( X ,0 = { A\ J \x,t)+A) J >{x-e u t 



U)u a i aU)i 



B\ j \x,t) = ± B\ J \x,t)+By>(x-ei,t) 



A I d(J)i 



(3.3) 
(3.4) 



Finally, a six-state basis was constructed from "transverse-smoothed" A- and B-fields that 
consisted of 

zero-gluon state: ... Qi(f) = l2, (3 5) 

one-gluon state: ... 62(f) = Aj (x + mj,t) , (3.6) 

,1=0 



4 



Gluon chains and the quark- antiquark potential 



Stefan Olejmk 



two-gluon states: ... 03(0 

... 04(0 
... 05(0 
... 06(0 

The antiquark is assumed to sit at Xq, the quark at Xn +Rej. Q\ is the zero-gluon operator, 02 
the simplest one-gluon operator, with one A-field put at different locations between the quark 
and antiquark, and 03_6 are simple two-gluon operators, containing two powers of the gauge 
field A. (In two-gluon operators 03 and 04 the interval of A-field insertions was extended to 
up to two lattice spacings outside the region defined by quark/antiquark positions.) 

Of course, one could use a larger set of more sophisticated operators and/or more variational pa- 
rameters, but we believe that the above choice allowed to fulfill the goals of the study in a clear-cut 
and convincing, even though only qualitative, way. 



R+1 R+1 o n 
= E J-Ay^xo + ne^OAV^Xo + n'e^O, 

n=—2 n'=n 
R+2 R+2 

= E L £a)%o + ^,04 (xo + «'e ; -,0, 

n=— 2 n'=n 
R-l 

n=Q 

= EL B i (xo + ^,012- (3.7) 

n=0#j 



4. Results 



The calculations outlined in Section g with the operators 0^ given by Eqs. (3.5-|3.7|), were 
carried out for SU(2) lattice gauge theory at coupling j8 = 2.2 on a 12 4 lattice volume, j8 = 2.3 
on a 16 4 lattice volume, and /3 = 2.4 on a 22 4 lattice volume. The operators 0^ depend implicitly 
on a variational parameter p , and matrix elements Tm and the gluon-chain potential V c h a in (ft) were 
computed for each R at twelve values of p, p„ = (n — l)Ap, 1 < n < 12, with Ap = 0.025 at j8 = 2.2, 
and Ap = 0.02 at j8 = 2.3,2.4. The choice of n which minimizes Vchain (ft) depends on both /3 and 
the quark separation R. For example, at j8 = 2.4 and R = 9, the optimal value was n = 8. All our 
data were always obtained from the optimal value of p for a given coupling and separation; below 
we will, with the exception of Fig. |5[ only report results for the largest coupling studied, j8 = 2.4. 



4.1 Gluon-chain potential 

In Fig. H| we compare the color-Coulomb potential with the gluon-chain potential. The latter 
was estimated using Eq. (2J5), while the former is given by Vc ou i = — log(7n), where T\\ is the 
zero-gluon to zero-gluon matrix element, independent of the variational parameter p. We also 
display the usual static quark-antiquark potential V trU e computed by standard methods from timelike 
Wilson loops with "fat" spacelike links. The data were fitted by the usual (constant + Liischer + 
linear term) function, the extracted string tensions were 0.158, 0.095, 0.069 for the color-Coulomb, 
gluon-chain, and true potentials, respectively. The inclusion of one- and two-gluon operators affects 
the string tension in the expected way: the Coulomb string tension, about 2.3 times larger than the 
true asymptotic string tension (at j8 = 2.4), goes down to the "chain" tension that differs from a asymp 
by 38% only. One can imagine that a modest improvement of our operator basis plus inclusion of a 
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few more constituent gluons would bring the string tension of the optimal variational gluon-chain 
state even closer to the true value. 



4.2 Effects of finite volume 

Figure [3] shows the dependence of the color-Coulomb potential and the gluon-chain potential 
on R for three different lattice volumes. While Vcoul bends away from linearity at R L/2, which 
is clearly a finite-size effect, V c hain seems completely insensitive to the size of the lattice. A natural 
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zero-gluon and one-gluon energies at (3=2.4 
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Figure 4: Energy of the zero-gluon and one-gluon states vs. quark-antiquark separation. 

interpretation of this result is that the long-range field does not exist or is greatly suppressed in the 
chain state relative to the color-dipole field of the zero-gluon state. 

4.3 Constituent gluon content of the gluon chain 

Finally we studied the constituent gluon content of our optimal variational gluon-chain state 
for different quark-antiquark separations R. In Figure |] we display the /^-dependence of the of the 
zero-gluon and one-gluon states at j8 = 2.4. The energies are given by — log(7n), and — log(r22) 




gluon content in the chain state 
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Figure 5: Zero-, one-, and two-gluon fraction of the norm of the variational state vs. quark-antiquark dis- 
tance R at j3 =2.2,2.3,2.4. 
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respectively. The Coulombic energy of the zero-gluon state, and the kinetic plus interaction energy 
of the one-gluon state become equal at R about 1 fm. 

One can estimate the gluon content also directly. If we write the minimal-energy variational 
state in our six-state basis as = Efc=i a k{R)\®k), then the zero-, one-, and two-gluon frac- 

tions are given by af, Ot^, and 1 — a\ — a%, respectively. The results are shown in Figure H The 
gluon content vs. R (expressed in physical units) turns out to be almost independent of coupling. 
The one-gluon content of the minimal energy state becomes equal to the zero-gluon content at 
about 1 fm, i.e. at the same distance at which the energies of zero- and one-gluon states equalize. 

5. Conclusions 

A simple variational calculation of a quark-antiquark state in a subspace of zero, one, and two 
constituent-gluon states yields its energy less than the energy computed from a zero-gluon state. 
This result is not surprising by itself - what appears nontrivial and was not guaranteed from the 
beginning is the following: 

1. The Coulombic energy of the zero-gluon state rises linearly with separation, albeit with string 
tension higher than the asymptotic tension of the QCD flux tube. 

2. The linearity of the potential survives addition of a small number of constituent gluons. 

3. A few constituent gluons tend to bring the string tension of the variational state down con- 
siderably, to a value closer to the asymptotic string tension of the QCD flux tube. 

4. One begins to see the formation of the gluon chain only at quark-antiquark separations of 
about 1 fermi. 
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